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Abstract 

The effects of quantum fluctuations on the magnetic field dependence of the spin-wave gaps 
in the lamellar antiferromagnet La2Cu04 are considered. Nonlinear corrections to the spin- wave 
spectrum are calculated to leading order in 1/S, where S is the localized spin. The nearest-neighbor 
exchange interactions between the Cu spins as well as the Dzyaloshinskii-Moriya interactions are 
taken into account. Using the experimental values of the components of the g-factor tensor, we 
get a satisfactory agreement with the experimental results for the field dependence of the gaps by 
Gozar et al. [Phys. Rev. Lett. 93, 027001 (2004)], and obtain consistent values of the in-plane and 
inter-plane coupling constants. The field dependence of the dispersion of spin waves propagating 
perpendicular to the Cu02 planes is also discussed. 

PACS numbers: 75.50.Ee, 75.30.Ds, 75.10.Jm, 74.72. Dn 
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I. INTRODUCTION 



Lanthanum cuprate, La 2 Cu04, the parent compound of the high-temperature super- 
conductors, is a layered weakly orthorhombic antiferromagnet (AF), with a Neel tempera- 
ture, T/v, of approximately 325 K.- The orthorhombic distortion is associated with a small 
tilt of oxygen octahedra around each copper ion, which also introduces the antisymmetric 
Dzyaloshinskii-Moriya (DM) superexchange interaction between neighboring spins. Because 
of the DM interaction, each Cu02 plane acquires a small ferromagnetic moment along the 
c-axis, perpendicular to the plane. The direction of the ferromagnetic moments alternates 
in adjacent planes, and, therefore, there is no net ferromagnetic moment in the crystal. 

The DM interaction, though small, has a strong impact on the magnetic properties of 
La2Cu04. It was established that a magnetic field perpendicular to the CuC>2 plane causes a 
first order weak- ferromagnetic transition (WFT).-^ The critical field, H c , of the transition 
depends on the DM coupling as well as on the interlayer exchange. In a non-stoichiometric 
sample (Tjy = 234 K), studied in Ref. [3], the critical field at low temperatures was H c « 4.8 
T. In Ref. the value H c m 6 T was obtained from Raman spectrum measurements for a 
crystal with T/v = 310 K, while in a recent paper- a significantly larger transition field of 
11.5 T was found from neutron diffraction studies for a sample with T/v = 316 K. 

Another manifestation of the DM interaction in La 2 Cu04 is the phase diagram in an 
in-plane magnetic field. The antiferomagnetic staggered moment in La2Cu04 is directed 
along the diagonal of the CuO plaquette (6-axis), the in-plane anisotropy being generated 
by the DM interaction. A theory^ based on the mean-field version of the Hamiltonian, which 
includes exchange and DM coupling between Cu spins, as well as an out-of-plane anisotropy 
term, predicts that in a field along the 6-axis one should observe two phase transitions: a 
spin-flop transition, at which the staggered moment jumps from the (b,c) plane into the 
(b,a) plane (a is perpendicular to b and c) at a field H, determined by the DM coupling 
energy, and a second transition at a higher field, when the staggered moment rotates into 



the c-direction. Such a phase diagram was indeed observed in Ref. 



Refs. 



7. Note, however, that in 



no spin-flop transition was observed. According to these references, the staggered 
moments rotate gradually from the b- to the c-axis. 

Gozar et ai— used Raman spectroscopy measurements to study the influence of a magnetic 
field on the in-plane spin-wave (SW) gap in La2Cu04. They showed that the gap abruptly 
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increases at the WFT, decreases with the increase of the field directed along the b- axis and 
increases with the field pointing along the a-axis. The observed field dependence of the gaps 
was described by classical theories^^ 1 ^ based on the model, which include anisotropic 
symmetric and antisymmetric intralayer coupling between the Cu spins as well as interlayer 
exchange. 

On the other hand, it is known that quantum corrections are very important when con- 
sidering the properties of layered cuprates, with spin 5 = 1/2.— Intensive calculations 
were performed for an isotropic Heisenberg two dimensional model with nearest neighbor 
interactions.— liiii^iiSiiZiiLiS Perturbative expansions in 1/5, series expansions, and Monte 
Carlo methods were used to go beyond the linear SW theory— iI&IS It appears that the 
renormalization factors obtained from the nonlinear SW (NLSW) theory to leading order in 
1/5 are sufficiently close to those which follow from more complicated calculations even for 
spin 5 = 1/2. The quantum effects renormalize significantly the mean value of the spin at 
T = 0, the SW stiffness, the in-plane and out-of-plane SW gaps at zero magnetic field, etc. 
In a recent paper,— the 1/5 NLSW theory was employed to calculate the effect of quantum 
fluctuations on the SW stiffness and other parameters of La2Cu04, in a model which takes 
into account the plaquette ring exchange. The renormalization factors which follow from 
these calculations and the fit of the experimental results by Coldea et al.,— appeared to 
be somewhat closer to unity than in the model with nearest-neighbor exchange only. In 
view of these results one should expect that quantum fluctuations are important also for the 
quantitative understanding of the SW spectrum behavior in a magnetic field. 

In this paper we use the NLSW theory to leading order in 1/5 to calculate the influence 
of quantum fluctuations on the magnetic field dependence of the in-plane and out-of-plane 
SW gaps in La2Cu04. We show that the quantum fluctuations alter the expression for the 
critical field of the WFT, renormalize all physical quantities, which enter in the equations 
for the gaps in a magnetic field, and change therefore significantly the values of the spin- 
spin couplings extracted from the experiment. The theoretical results are in satisfactory 
agreement with the experimental findings. The theory also predicts a non-trivial effect of 
the magnetic field on the SW propagating in the c-direction. 

The paper is organized as follows. In section II we present the model of magnetic inter- 
actions in La2Cu04 used in our calculations. In section III we derive the field dependance 
of the gaps for the field in the c-direction in both the linear and the nonlinear SW approxi- 
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mations. We also consider the field dependence of the dispersion of SW, which propagate in 
the c-direction. In section IV the NLSW theory for the SW gaps in a magnetic field along 
the b- and a-directions is developed. In section V we compare our NLSW results with the 
experiments by Gozar et a/.- Finally in the Appendix the Green's function are used to derive 
a general expression for the eigenvalue matrix of a model described by the Hamiltonian of 
the type (I14p . The Appendix also presents some useful relations between the averages of 
the magnon operators. 

II. THE MODEL 

Our Hamiltonian is written as^ 

n = j [s* • s,- - as?s*] + J2 D ^ ■ S * x S j 

<ij> <ij> 

+ J ± J2 S* • S* - X> fl H ■ g • (1) 

<ik> i 

Here the first term describes the anisotropic exchange coupling between nearest-neighbor 
spins in the Cu02 layers, the easy plane being the crystallographic (a, b) plane. The 
anisotropy factor a in L^CuO^ (LCO) is small, of order 10~ 4 .— The second term is the 
DM interaction, which leads to the canting of the spins out of the plane. The third term is 
the small interlayer exchange coupling. The sum is over the (two) nearest neighbors along 
the crystallographic c-axis, which we choose as the z-axis. Finally, the last term in Eq. (jTJ) 
describes the Zeeman energy in an external field H, g being the anisotropic ^-factor; fiB is 
the Bohr magneton. 

The DM vector Dy is perpendicular to the Cu-Cu bonds and changes sign from one bond 
to the next one. In the orthorhombic crystallographic axes (a,b,c), where a and b point 
along the diagonals of the Cu plaquette, the vector D^- can be written as 

B u+Sl = (A AO) 

B u+S2 = (0,-0,0), (2) 

where 5\ and 5 2 are directed along the two orthogonal bonds of the plaquette. Using these 
DM vectors, the DM term in the Hamiltonian transforms into 

Djj ■ Sj x Sj = D ■ Sj x Sj 

<ij> <ij> 
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FIG. 1: The four copper spins A,B,C,D in a unit cell. Ordered spin directions are shown by arrows. 
The horizontal lines refer to the CuO"2 layers. The global reference axes z and y point along the 
orthorhombic crystallographic axes c and b respectively. 



with the vector D of length D pointing along a. The axes (x, y, z) are along (a, b, c). 

The second term in the r.h.s. of this expression does not contribute to the gap. It only 
changes slightly, to order D/J, the SW dispersion.— We will therefore neglect it in what 
follows. Thus, for our purposes the DM interaction can be described by the vector D, 
directed along the a- axis in each plaquette. 

The magnetic unit cell in the above model contains four spins, A, B, C, D, in two inequiv- 
alent layers, Fig. 1. Due to the DM interaction, these spins cant into the c-direction, with 
canting angles, with /i = a,b,c,d for the spins in the unit cell. These angles depend on 
the direction and on the strength of the magnetic field. In zero field, they satisfy the relation 
if) a = ipf) = —ip c = —ipd-, so that the magnetic moment of the unit cell is equal to zero. To 
calculate the SW spectrum, it is convenient to rotate the spins from the above global crys- 
tallographic axes (x,y,z) to local axes (£, ?7, C), with ( along the sublattice magnetization. 
When the magnetic field is in the 6-direction, and smaller than the spin-flop field, or in the 
c- direction, the spins rotate in the (y, z) plane. Then the transformation equations from 
the crystallographic axes to the local ones for spins a and b in the sublattices A and B are 
(see Fig. 1) 



I r~) \ \ " / HI QZ QZ QX \ \ , / HI QZ QZ QX \ I o\ 




qy 



S2 sin ip a + S% cos i/) a , S z a = -S2 cos tp a + S% sin tp, 
S% sin - Si cos if) b , S b = Sjj cos ip h + Sf sin ip b . 



(4) 



Analogous equations hold for the spins c and d in the second layer. 
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III. FIELD ALONG THE c-AXIS 



The canting angles for this direction of the field are the same for each pair of the spins 
in a plane, i.e. ip a — ip b = ipi, and ip c — ip d = ip 2 . ipx and ip 2 , and, hence, the magnetic 
moments of the planes differ if H is smaller than a critical field H c , which in the classical 
approximation is equal to^ 



UJ±S 



(5) 



At H = H c the weak-ferromagnetic transition (WFT) happens, and the magnetic moments 
of all planes became the same.- At fields of order H c , the angles ip\ and ip2 are small, of 
order 10 -2 . 

In the local framework, the Hamiltonian (CQ) has the form 



7i = Tii + H.2 + TC±. 



(6) 



Here 



Hi =1 , 2 = -(Jcos2^ + Asin2^ + «Jsin 2 ^) £ S^^S^j) 

<ij> 

+ I(jcos 2 ^ + ^sin2^-^) + 

A z 2 Kij> 

+ I(j s in 2 ^-^sin2^ + ^) £ [ S+ W^"0') + S~{i)S + U) 



- g c fi B H sin [J2 S<(i) + £ S^j)] + g cf i B H cos ^ ^(i) - £ S ' n U) 
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aJ 



+ (Jsin2^ - D t cos 2^ - — sin 2^) £ ['S'ty'W) " ^(O^O 



where .D/ 



-1VD, 



(7) 



(8) 



and labels the spins in the sublattice A (B) 
The interlayer Hamiltonian transforms into 

1 r 



n± = -j± "£{s<(i)s({k) cosVi,2 



<ik> 



S + (i)S + (k) + S-(i)S-(k) 



cos 



2 ^1,2 



1 r 



S + (i)S-(k) + S-(i)S + (k) 



■ 2 ^1,2 . 
sm — r- + 



S c (^(£;) - S^)S ? (£;)] sin^ 2 }, (9) 



where ^1,2 = V'i + V*2- We next make the Dyson-Maleev transformation to the bosons 
operators a, b, c, d for the four spin operators in a unit cell. In the local framework, this 



transformation is the same for all spins. It can be written as 



S+ = V2Sa„ S~ = V2Sal(l 



2S 

S C , = S- afa, (10) 
where the operators a M are defined as a.\ = a, a 2 = b, a% = c, = d. 



A. Linear spin-wave theory 

It follows from Eqs ©, O, @ that the classical ground-state energy per spin, E , is 

E = ^[-JS 2 z{cos2ip 1 + cos 2ip 2 ) + DS 2 z(sm2ip 1 - sin2^ 2 ) 

— 2g c fi B HS(sinipi + sin ip 2 ) — 4Jj_S' 2 cos(V>i + ^2) — aJS 2 z(sm 2 ipi + sin 2 ^2)], (H) 

where z = 4 is the number of nearest neighbors in the plane. The minimization conditions 
for the energy E Q with respect to ipi an d 4>2 yield 

JSz{2 — a) sin 2ipi + 4Jj_5 , sin(^i + ip 2 ) + 2DSz cos 2ipi — 2g c \x B H cos ipi = 0, 
JSz(2 -a) sin 2ip 2 + 4 J± S sin(V>i + ip 2 ) -2DSz cos2il; 2 - 2g c /i B H cosip 2 = 0. (12) 

The solution of these equations gives ipi and ip 2 at fields H < H c . Given that ip± and ip 2 are 
small, and neglecting terms of order a and J z /J, we get 

g c li B H - DSz g c jj, B H + DSz 
^ = 2JSz ' ^ = 2JSz • (13) 

With these angles, the terms in Eq. ([7]) which are proportional to S v , and are therefore 
linear in the boson operators, cancel out. 

The bilinear spin-wave Hamiltonian, which follows after the transformation ( flOj) is per- 
formed, can be written as 

H {2) = £{V(q)aI(qH(q) + -B^(q)[al (q)at (-q) + a M (q)a„(-q)]}, (14) 
q A 

where the matrices A(q) and B(q) are given by 

A u = A 22 = JSz{\ - 2R X ) + g^H^ + 2J ± S, 
A 33 = A 44 = JSz(l - 2R 2 ) + gcf jL B Hifi2 + 2J L S, 
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Al2 


= A 


A^ 


= A 


A u 


= A 




= B 


-B34 


= B 




= B 


B n 


= B 



JSzjq^Rl + — j, 

2:-; — A13 = A14 = 0, 



= JSz^lA - R 2 - -J, 
2:; = 2J±Scz, 

22 = -B33 = i? 44 = Bi 3 = B 3 i = B 2 a = -B42 = 0, (15) 



while 



R ^ = H^-T)= A{JSzf ' (16) 
cos q x + cos q y 

7q = g ' c ^ = cos ^- ( 17 ) 

The wave-vector components in the plane and along the z-axis are measured in units of the 
corresponding intersite distances. In Eq. ( fl5|) we neglect corrections of order J±ip? and aipf 

(1=1,2) 

see also the Appendix of this paper) that the squares of the 



and 
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It was shown in Ref. 
spin- wave energy, u 2 , are the eigenvalues of the matrix 

M = [A(q) + B(q)] X [A(q) - B(q)]. (18) 

Since the magnetic field is small, g c /j,BH <C J, it is not expected to influence essentially the 
stiffness of the spin-waves propagating in the (x, y) plane. In contrast, the effect of the field 
on the spin waves propagating in the z-direction may be strong. We consider, therefore, 
only the dispersion of the spin waves propagating in this direction, and put 7 q = 1. The 
spin- wave spectrum at H < H c , which follows then from Eqs. (fT5]l and ( TlBl) . is given by 



uf n = 5 2 + w 2 ±^(g cf i B H6) 2 + w*cl 

u 2 out = A 2 + (g c fi B H) 2 + w 2 ± J(g cf i B H5) 2 + w*c 2 z . (19) 



Here Ui n and u out are the frequencies of the in-plane and out-of-plane spin waves, respectively. 
The zero-field in-plane, 5, and out-of-plane, A, gaps are: 



5 = DSz, A = JSzV2a, (20) 
8 



and w 2 = 4zJJ±S 2 . The lower (upper) sign in Eqs. (TT9]) gives the frequencies of the in-phase 
or acoustical (out-of-phase or optical) oscillations in adjacent planes. 

At H = H c , the field which favors alignment of magnetic moments in adjacent planes 
overcomes the interlayer coupling which favors alternation of the moments. Then the stag- 
gered moments of all the planes align in the same direction, accompanied by a 180° rotation 
of the spins in half of the planes. Hence, the canting angles at H > H c are 

g c fi B H + DSz 

*> = *-^= 2JSz ■ (21) 

The difference AE , between the ground state energy E at H < H c and that at H > H c , is 

= g c » B HDS _ 2 _ 
° 2J v ; 

The expression ([5]) for H c follows from the relation AEq = 0. 

Since the two adjacent planes are equivalent, only the two modes which describe the 

in-phase oscillations in the planes are relevant. One gets for these modes 

u in = + 9c^bH5 + w 2 (l - c z ), 

uo 2 out = A 2 + (g c fi B H) 2 + g c fi B H5 + w 2 (l - c z ). (23) 
The spin- wave gaps which follow from Eqs. ffl9|) and (123]) at c z — 1, coincide with those 



obtained in different ways in Refs. |9h12 



B. Nonlinear spin waves 

We now consider the effect of terms which are cubic and quartic in the boson operators. 
It follows from Eqs. (J7j) and (Q that all the terms which are linear in S v , except those which 
are proportional to H, cancel out owing to Eq. ffT2]) . Thus, the only terms odd (cubic) in 
boson operators that remain in the Hamiltonian are proportional to H. For one plane and 
for small fields g^sH -C JSz we have 

= g c ^ B H\U- 1 - J2 WMh - &}) + &}&i(oJ - *)]■ (24) 

v Zb z <ij> 

The fourth-order terms of order 1/S come from the exchange and DM interactions, and from 
the out-of-plane anisotropy. They are 

Hf = -J{1 - 2^ 2 ) £ a\ ai b)b 3 - i{\ - ^) 53 a\b](a\a t + b}bj) 

<ij> <ij> 
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53 a l a l«i & i> 

<*i> 



^(4) 

= 2 



[4ajaj6]6j + a|a|aj(6] — &,) + (a| — aj)&j&j&j], 

<ij> 

Ja 



53 [(a* - aJ)&J-6j-6j + (6j - b))a\a\ 



O; 



(25) 

(26) 
(27) 



<r/> 



The 1/S* corrections to the spin wave spectrum come from the first-order contribution of 
Ti^ and the second-order contribution of Ti.^. The last contribution is of order H 2 , and it 
renormalizes the H 2 terms in the out-of plane spin gap, Eqs. (TT9l) . ( 1231) . It was calculated 
in Ref. 25. 



To treat the quartic perturbation, we truncate all four-operator terms by contracting out 
pairs of operators in all possible ways and discarding the non-Hermitian terms.— One finally 
arrives at an effective bilinear Hamiltonian, of the form shown in Eq. (JHj), with the matrix 
elements A^ v and B^ v depending on the following averages of the magnon operators: 



v 



< aUi >= 53 «i«q ; f = 51 < fl q & -q > 7q> 



V = 51 < a q fe q > 7q> ^ = 53 < a Q a -q > • 
q q 

The matrix elements A^ u and B^ u , with \i — 1, 2; v — 1, 2 can be written as 

^11 = ^22 = JSZZ, 



(28) 



1 - 2i? 1 (l + -L) + ||] + #^ + 2J 2 SZ m , 



Ala 
#11 



r~ 1 a 

A 2 i = J5z7 q [ J RiZ m - —(A + 2r]) + Z m - 
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#21 = ^7q^c 1-^1 1 



2( 1 + l 



Here is the tilting angle renormalized by quantum fluctuations (see below), 



D 
1 



(29) 



(30) 



Z c = 1 — u/S — ^/Sis the renormalization of the SW velocity,— 1 ^ and Z m = 1 — v/S is the 
renormalization factor of the average spin. 

The equations for the matrix elements A^ v and B^, with /x = 3, 4; i/ = 3, 4 follow from 
the above equations by substituting ip 2 for ip 1 . The mixed elements B u and -B 23 are 



Bi4 — B 2 3 — 2J±SZ m c z 



(31) 
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Note that averages like < OjCj > appear in the course of the renormalization. They are 
small, since i and j belong to different planes. When deriving Eqs. fl29l) we neglected also 
terms of order of r]if) n or \ip n , since (as shown in the Appendix) A and r\ are of order of ipf. 
The renormalization factors £ and v are known^^ to be equal to 

lw 1 ^ „ 1n „ , U 7, 2 



nEh=^ - 1 = 0-197, C = -,E ~7== = -0-276- (32) 



This gives Z c = 1 + 0.158/25 1 = 1.158. More accurate values of Z c , obtained from expanding 
to order 1/S 2 or from Monte Carlo simulations, are also available.— 1 ^ Nevertheless, for 
consistency we use in what follows this value of Z c . The factors 77 and A are given in the 
Appendix. It appears that the spin- wave frequencies depend only on the sum rj + A, which 
according to Eqs. (I29I) and (1A.10I) is equal to 

A + ?7 = -(# 1 + -)£. (33) 

To calculate the renormalized tilting angle we add to the classical ground-state energy 
( TlTl) the (1/S) correction, which follows from the bilinear Hamiltonian f THl) after averaging 
the products of the boson operators. This yields for the quantum corrected ground-state 
energy, E: 

E = ^{ (l - - ) [- JS 2 z(cos 2^i + cos 2^2) + DS 2 z(sm 2^i - sin 2^2)] 

- 2g cf i B HS(l - |)(sinVi + sin^) - 4J ± S 2 (l - ^) cos(</>i + tfa)}. (34) 

We neglected here terms of order rjipi, as well as small corrections which come from the 
out-of-plane anisotropy. 

Minimization of the energy ( 1341) gives 

^ = 2JZ^ ± 2J' (35) 



in agreement with previous calculations for an AF with .D = 0.— >2£ Note that the contri- 
bution to the canting angle caused by the DM interaction is the same as in the classical 
limit. 

When calculating the in-plane spin- wave energy from Eqs. ({TBI and (|29|) . one finds a 
contribution to the gap which is of order 1/S" and proportional to H 2 . This contribution 
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should be cancelled out by the contribution of the cubic terms in the Hamiltonian, since the 
H 2 term in the in-plane gap is forbidden by symmetry. Discarding this term, we get 



col = s 2 + w 2 ± yW#) 2 (i - fp 2 + 



u 2 mt = A 2 + Z c {g c ^ B H) 2 + w 2 ± ] j {g^ B H) 2 (l - ^f~5 2 + w*c 2 . (36) 

Here 5 = Z m 5 = (1 - v/S)5, A = Z m A, w = Z w w, where Z vl = l-v/3- £/2S « 0.88. As 
noticed above, the renormalization of the H 2 term in the out-of-plane wave was calculated 
in Ref. 25. 



We see that the quantum fluctuations effectively increase the ^-factor in the c-direction. 
This is in contrast to the result,— obtained in the limit of large N, N being the number of 
the components of the spin. 

Above the WFT, when the staggered moments in all planes align in the same direction, 
the angle ip2 is given by Eq. (|35|) . while ipx is given by ipi = n — ip 2 - It follows from Eq. (|34|) 
that the difference, AE in the ground state energy before and after the transition is 

^=2=^(1-^-2^(1-1). (37) 

The condition AE = then yields the transition field, 

H e = UZmJ ' S = g(l + |). (38) 



The above expression for H c follows from the following simple qualitative arguments. At 
the WFT the magnetic moment increases abruptly by the value < S z > D / J = Z m SD/J 
[see Eq. ( 1351) ]. and the gain in the magnetic energy per spin is < S z > g^sH D / J . The loss 
in the exchange energy is 2J Z < S z > 2 . Equating these energies one arrives at Eq. (|38|) . 
The SW spectrum at H > H c is 



= d 2 • Ud(\ • ir 2 {\ -r ;: ). 



s 



oo 2 out = A 2 + Z C H 2 + H~8(l - |) + w 2 (l - c z ). (39) 
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IV. FIELD IN THE (a, b) PLANE 



A. Field along the 6-axis 

When a magnetic field is applied in the direction of the staggered magnetization, i.e. 
along the fe-axis, an unusual spin-flop transition happens at the field H = 5.— At this field 
the staggered moments rotate from the (b, c) plane to the (a, c) one, forming an angle 
with the c-axis, which decreases with the increase of the field. At a higher field, equal to 
H = (A 2 + 2w 2 — S 2 )/S, the staggered magnetization points in the c-direction.- 

We consider in this paper only fields smaller than S, when the spins lie in the (b, c) plane, 
which we defined as the (y, z) plane. As in the case of a field along the c-axis, the DM 
interaction and the field cause the canting of the spins out of the (a, b) plane. The canting 
angles of the spins a and b, however, are no longer equal to each other. 

After the rotation to the local framework is performed, we obtain the Hamiltonian for 
the first plane (Fig. 1) as follows: 

Hi = H ev + 7~Lodd- (40) 

Here 

H ev = -JcosW* E ^)^(j) + ^cos 2 ^ £ [s + ws + (j) + s-ws-(j) 
<ij> 1 1 <ij> 

+ |sin 2 ^ J2 [5 + «5-(j) + 5-(z)5 + (j)] + gwBHifxxfr- cos^E 5 ^) 

1 1 <ij> i 

+ Dsin^ E [S<(i)SHj)-\(S + (i)-S-(i))(S + (j)-S-(j))~ 

<ij> 

- jEN-^lN-^')). ( 41 ) 



and 



H oM = (JsinVu + ^cosVwO E [S c (i)S'U)-S'(i)S<Vj 

<ij> 

+ g b /i B H (sin ip a - sin if) b ) ^ S v {i) 

i 

U [S C (i)S"U) Bin^aCOsVft - S^ijS^j) sin^cos^a], (42) 



<ij> 



where ip a b = ifta + The Hamiltonian for the second plane can be obtained from Eqs. ( )4T]) 
and (T42l) by replacing ip a and ip b by ip c and ipd, and changing the sign of D. Finally, the 
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interlayer coupling is 

H± = J± Y / {sH^{k)-[S<(i)S^k) + S r 'ii)S r '(k)]cos('i/; a + ' l p d )} 

<ik> 

+ j± E {sHj)sHk) - [sHj)s c (k) + s"0W*)] «»(^ + v c )}- 

<jfc> 

It follows from the above equations that the classical ground-state energy is 

E = ^{-JS 2 z[cos(ip a + ip b ) + cos(ip c + ip d ) 

+ a(sin ip a sin %p b + sin ip c sin ip d )] - DS 2 z[sin(ip a + ip b ) - sm(ip c + ip d )\ 

+ g b fisHS (cos ip a — cos ip b + cos ip c — cos Ipd) 

- 2J ± S 2 [cos(ip a + ip d ) + cosOfc + ipc)}}- 

The minimization conditions with respect to ip^ yield 

DSz(gbf^BH — iJ±S — aJSz) 



Ipa 
Ipb 



-Ipc 

-4>d 



A 2 + 2w 2 - (g bt M B Hy 
DSz(gb/J>BH + iJ±S + aJSz) 



(43) 



(44) 



(45) 



A 2 + 2 W 2 - tg bf JL B H)* ' 
When deriving Eqs. ( )45l) we assumed that ip^ is small, i.e. we consider magnetic fields that 
satisfy the inequality g b [iBH <C A. 



1. Spin gaps in the one-layer model 

As we show below, the interlayer coupling almost does not affect the spin-wave gaps. 
Hence, when calculating the gaps, one may put J± to zero, and consider spin waves only in 
one plane. Performing the Dyson- Maleev transformation of the corresponding terms in the 
Hamiltonian (14TI) . one gets the bilinear and quartic terms in the Hamiltonian as follows: 



n {2) = Y.[Js{i-^)-Ds^p ab 



+ E [^(1 - - ^ - ^f] (a ibj + ajftt) 



DS ~ 



aJS 



<ij> 



1 JS 

+ o E \~^~^lb + DSip a b + aJS (aib] + a\bj) + gbfi B Hj2( a i a i ~ (46) 



2 ^ L 2 



ww = -^E[( 1 -%4^ + k 1 + % a H(^+^) 
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ib 2 

- ^(atytybj + a\a\aibj) 

+ Yl{ D $ab a i a ibjbj - -(D^ + aJ)[(ai - a\)b]b]bj + ajala^bj - &])]}. (47) 

u 

Here ip ab = ipa + 4>b is the renormalized tilting angle. To calculate it one should add, as 
in Section III, to the classical ground state energy (j44j) the quantum correction, i.e. the 
averaged SW energy, and then minimize the full energy. 

The terms linear in the boson operators in Eq. ( I4"2l cancel out when the angles given 
by Eqs. (f45l) are used. The cubic terms partly cancel out as in the case when H is along 
the c-axis. The remaining cubic terms are of second order in H and in the other small 
parameters of the problem, and are therefore irrelevant. Hence, the quantum corrections 
can be calculated in the Hartree approximation as was explained in the previous Section. 
This gives an effective bilinear Hamiltonian, with the coefficients given by 

£ \ a 



A n (H) 
A 22 (H) 
A 12 

fin 

B12 



JSzZ r 



/ C \ Ut 1 



An(-H), 
JSz 



2 

Jz 

"T 



R - -(X + 2r]) + a 



R 



1 



1 



Here 



R = ^ab(— + -r), 



s 



1 + 



s 



2 ' J 

and the angle ipab renormalized by quantum fluctuations is 

2aJD(Z c Sz) 2 



ab 



(Z C A) 2 - (g bf M B H) 2 
Equations ([50]) and ([39]) yield 

a(Z c 5) 2 [Z 2 A 2 -2(g b ^ B H) 2 } 



R 



+ 0(H 4 



(48) 



(49) 



(50) 



(51) 



[Z 2 A 2 - (g b ^ B H) 2 ] 2 2J 2 
Thus, the field dependence of the gaps comes only from the explicit dependence on the field 
of the matrix elements An and A 22 . We get for the in-plane, VL ini and out-of-plane, Q ou t, 
gaps the following expressions 
lr.o 



n 2 



tt 2 



5 + A + 2(g b u, B H) 2 - V(A 2 - 5 2 ) 2 + 8(g bf i B H) 2 (5 2 + A' 



5 + A + 2(g b ^ B H) 2 + V(A 2 - 6 2 ) 2 + 8(g bf i B H) 2 (S 2 + A 2 



(52) 
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The in-plane frequency Qf n changes sign at the field gbfiBH = 5, when the spin-flop transition 
takes place. 

2. Effect of interlayer coupling 

The generalization to a three-dimensional crystal is straightforward. Even though ip^ 
now depends on the interlayer coupling w, the quantity R is, as before, equal to D 2 /2J 2 
up to terms of order H A . The in-phase gaps therefore do not contain w, and as before they 
are given by Eqs. (1521) . The spectrum for spin- waves propagating along the z-axis in small 
fields Qbl^BH C A is given by 

2/ \ r>2 , ~2|"i 4(# 6 /^ B if) 2 -i 

^in(^) = Q in + W [1- —^—~ r \(l±C z ), 



"IM = rf out + w 2 \l + 4( |fff 2 )S ] (1 ± c z ). (53) 



Here the upper (lower) sign corresponds to the in-phase (out-of-phase) modes. 
B. Field along the a-axis 

In this case the field and the DM coupling do not interfere. Hence, the spectrum is the 
same as one would have for an AF in a field perpendicular to the easy anisotropy axis. 25 
The only frequencies relevant for the acoustic spin waves in this case are given by 

uUVz) = $ 2 + Z c (g af i B H) 2 + w 2 (l - c z ), 

"IM = A 2 + w 2 (l-c z ). (54) 

V. DISCUSSION: THEORY VERSUS EXPERIMENT 

Consider first the field in the c-direction. The value of g c in La2Cu04 is not well deter- 
mined yet (see below). We therefore exclude g c from the expressions for the in-plane gaps, 
Eqs. (EED and (ESJ), using Eq. (EHJ), and obtain: 



~5 2 -w 2 (sjH 2 /H 2 + c 2 z -l), H<H c , 
^1 = & + w 2 (H/H c ) + w 2 (l - c z ), H > H c . (55) 
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These expressions differ from their classical counterparts only by the parameters 5 and w, 
which should replace their classical values. Thus with the proper choice of the parameter w 



(S and H c are measured in the experiment 
a/.- in the same way as was done in Refs. |5j, 
duf n (H > H c )/dH is given by 

dul(H > H c ) _ 
dH 



one can fit the experimental data by Gozar et 
12 . According to Ref . {], the derivative 



11 



and 



-Q35 ( meV f 

H r T 



(56) 



There is an uncertainty in the value of H c , because of the hysteresis observed in the first-order 
WFT. We choose the value H c = 6.3 T, at which the gap ceases to decrease with the increase 
of the field.- With this value of H c we obtain from Eq. (1561) that w 2 = AZ^JJ ± = 2.19 
(meV) 2 . Given that J = 135 meV,— we get: 

J ± = 5.3 ■ KT 3 meV, a± = J±/J = 3.9 ■ 1(T 5 . (57) 

The ratio w 2 /H c can also be obtained from the relation 

w 2 /H c =~5g c ii B {l-i/S), (58) 

if the g-tensor is known. 

According to the calculations performed in Ref. 
are: g c = 2.45, g a = 9b = 2.11. Similar values were later obtained in Ref. 
experiment gave somewhat lower values: g c = 2.3, g a = 9b = 2.08,— or g c = 2.24, g a = gb = 
2. 06. 31 In what follows we use the last values of the g-tensor. Given that 5 is equal to 2.16 
meV,- we have: w 2 /H c = 0.43 (meV) 2 /T. This value is somewhat larger than that, which 
follows from the fit of the field dependence of the gap, Eq. ( |56l) . The discrepancy might be 
caused by the insufficient accuracy of the 1/S expansion, since the factor 1 — £/S — 1.55 is 
significantly larger than one. 

With 5 = 2.16 meV (2.3 meV according to Ref. |23| ) we obtain the DM coupling as: 
D = S/2Z m = 1.78 meV (1.90 meV). This value is by a factor two larger than that obtained 



3dL the principal values of the q-tensor 



12. The 



in Ref. U Note that the renormalization of both the in-plane and out-of-plane gaps is given 



23l . For more discussion on this subject see Ref. 



by Z m rather than Z c used in Refs. |5| and 
29. 

The dispersion of SW propagating in the c-direction also depends strongly on the field. 
We have at H < H r 



"1(<U) - ^n(0) = w 2 (^l + H 2 /H 2 - ^cl + E 2 IE\ 



(59) 
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The dispersion decreases with the increase of the field, and at the critical field, H = H C} it 
jumps to the value in zero field, 

^ n (q z )-u J 2 n (0)=w 2 (l-c z ). (60) 

At fields larger than H c the dispersion does not depend on the field. 

When the field is in the 6-direction, the measured in-plane gap decreases with the increase 
of the field as^ u 2 n = 5 2 — jbH 2 , with 7^ = 0.025 (meV/T) 2 . This dependence follows from 
Eq. (152!) for fields small in comparison with the out-of-plane gap A. In this case we have 

1 + 3(5 /A) 2 . 2 

lb = ~ -/ (gbVB) 1 • (61) 

1 - {5/ A) 2 V ' K ' 

Given g b = 2.06, 5 = 2.16 meV, and the above value of 75, we get A = 5.3 meV, in good 
agreement with the neutron measurement result A = 5.5 meV.— It follows from this result 
that the out-of-plane anisotropy parameter a is 

a = A 2 /8J 2 Z 2 m = 5.2 ■ 10" 4 . (62) 

The value of a is by an order of magnitude larger than a±. This implies that T/v in lanthanum 
cuprate is determined mainly by the out-of-plane anisotropy. 

Finally, we consider the gap when the magnetic field points in the a-direction. The 
experimental findings were fitted in Ref. ^ by the relation u 2 n = 5 2 + "y a H 2 , with 7 a = 0.015 
(meV/T) 2 . The theoretical value, which follows from Eq. (|5^|) . with g a = 2.06, is somewhat 
larger: 7 a = 0.0165 (meV/T) 2 . Note that in the model, which includes ring exchange, Z c 
was found to be Z c = 0.96M With this value of Z c one gets 7 a = 0.014 (meV/T) 2 . 

In conclusion, we calculated the effect of quantum fluctuations on the SW spectrum in 
La2Cu04 in an external magnetic field. With the experimental values of the g-tensor, the 
in-plane gap, and the WFT critical field, a satisfactory agreement of the theory with the 
experimental findings was obtained. Given the renormalized values of the parameters, which 
determine the field dependence of the SW gaps, we got new values of the in-plane DM and 
inter-plane couplings, and of the in-plane anisotropy a. 
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APPENDIX 

Define the Green's functions 

GW(q, *) =« ^(q, *), «l(q) »> F^(<h t) =« aj(q, t)at(-q) » . (A.1) 
These functions obey the following equations of motion: 

i^j-G^fat) = 6(t)S lu/ + « [a M ,?<],c4(q) >> 

ijF^t) = « [ai(-q),^],4(q) » . (A.2) 
Given the Hamiltonian (j!4p . one gets for the Fourier transform of the functions G and F 

(wI-A(q))G(u;,q)-B(q)F(a;,q) = I, 

(cuI + A(q))F(^q) + B(q)G(^q) = 0, (A.3) 

where I is a unit matrix. 

The equations of the eigenvalue problem can be written as 

{uI-A)\G > -B\F > = 0, 

(ojI + A)\F > +B\G > = 0, (A.4) 

where \G > and \F > are columns of the corresponding matrices. 
This gives 

u\G-F> = (A + B)|G + F>, 
u\G + F> = (A-B)\G - F > . 



Hence, 



(A.5) 



u 2 l- (A + B)(A-B) =0. (A.6) 
Thus, the squares of the spin- wave energy, u> 2 , are the eigenvalues of the matrix— 

M = [A(q) + B(q)] X [A(q) - B(q)]. (A.7) 
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The quantities rj and A, denned in Eq. ([28]) . are related to the Green's functions as 

n J - 

- I dw^F 11 (u,q). (A.8) 

7T J 



V 

7T 



A = 

7T 



Together with Eqs. flA.3j) . this leads to the relations 

A 12 _ 7* 



V 



2J^(l- 7 2) 3 / 2 ' 

Al2 7o 



Thus, we have 

77 -I- A = 

JSz 



V + \ = 4^Z- (A.10) 
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